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Abstract 

The hadroproduction of D s mesons is discussed within 
the framework of the modified Quark-Gluon String model 
(QGSM) taking into account the decays of corresponding 
S-wave resonances. A description of the existing experi¬ 
mental data on inclusive spectra of D s - mesons produc¬ 
tion in irp - collision and asymmetry of leading (Dj) and 
nonleading (Df ) mesons in £ _ p- interaction is obtained. 
The predictions for cross sections and leading/nonleading 
asymmetry in irp - and T,~p- collision for different initial 
momenta are also given. 
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Hadroproduction of charmed particles is now being inves¬ 
tigated both experimentally and theoretically. On the one hand 
the methods, available for light hadrons such as quark models, 
symmetries and so on, can be applied to hadrons, containing the 
charmed quark. On the other hand, the presence of c- quark inside 
the charmed hadrons allows one to use QCD methods. 

There is wide variety of theoretical models which are more 
or less satisfactory applied to description of charmed particle 
hadroproduction (see [|lj and references therein). 

Among these models the Quark-Gluon String model (QGSM) 
was successfully used for the description of many features of mul¬ 
tiparticle production in hadron-hadron collisions both for light 
and charmed hadrons. 

The modification of the QGSM taking into account the con¬ 
tributions from decays of corresponding S'-wave resonances was 
developed in [|], 

In the previous papers 0.0 this approach was used for the 
calculation of charmed hadron (meson and baryon) spectra, pro¬ 
duced by different hadron beams (tt, p, £ and S) without taking 
into account the charmed sea contribution. 

The charmed sea contribution in QGSM for the D mesons 
hadroproduction were considered in [f|, |J to describe the experi¬ 
mental data of leading/nonleading D— mesons asymmetry. 

The predictions for inclusive spectra of D s - mesons produced 
by different hadron beams in the various modifications of QGSM 
without taking into account the contribution of charmed sea were 
calculated earlier in il- However in these papers, due to the 
absence of the measurements, there was not made the comparison 
of the model calculations with experimental data. 

In this note we analyze for the first time the recently ap¬ 
peared experimental data 0 - [|1^] on inclusive spectra and pro¬ 
duction asymmetry of the leading/nonleading charmed-strange 
Dj(sc) - and D+(cs)~ pseudoscalar mesons. The analysis has 
been performed in the framework of modified QGSM ||2| taking 


1 


into account the contribution of vector D* resonance decay into 


D s through D* —> Ah7 mode 


The formulae for inclusive spectra of D* mesons hadroproduc- 
tion for different beams and the expressions for quark and diquark 
distribution functions in the initial hadrons were given in ||, [4j . 
Here we used the parametrization of charmed sea in and 
given in Q. 

In the present paper we use the following parametrization of 
the quark and diquark fragmentation functions into D s 


mesons: 


G?(z) 


= (1 


NA-a v ,(0)+2(l-a H (0))+2(a fl (0)-a 9 ,(0)) 


( 1 ) 


G?(z) = Gi-(z) = G?{z) = G?(z) = Gf{z) = 

Gf(z) = G^(z) = G^(z) (2) 

G?‘ 7 (z) = (1 - z) x ~ a ^ 0 \l + a f z 2 ) (3) 

(z) = (1 - z) x - ^(o)+2(i-Mo)) ( 4 ) 

G° ( 7 J Dt \ Z ) = 1 _ ^)A-a v (0)_ ( 5 ) 

a 0 s 

Gdl ( z ) = (1 — z^~ a ^ (0)+ 2 ( 1— (0))+(oh(0)—(0)) (Q) 

G£ (z) = G d J (z) = G% (z) = G$ (z) (7) 


(') = (1 ~^+<*r( 0)—2ajv(0)+ctii(0)—^( 


1 + a,\z 2 (1 — zf 


( 8 ) 


qD+ ^ ^ ^ A—ay,(0)+2(l—ajv(0))+2 (or( 0)—a v (0)) 


(9) 
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where A = 0.5, or( 0) = 0.5, «7v(0) = —0.5. Under these 
calculations we consider two values of the intercept of charmed 
cc- trajectory a^( 0) = —2.18 and a^( 0) = 0.. 

Now we turn to the comparison of the existing experimental 
data 0- 0 for Z) s -mesons produced in n p and E p collisions 
with the QGSM calculations in the framework of the model under 
consideration using fragmentation functions (0) - (|^). The values 
of free parameters were obtained by fitting to the experimental 
data. 

In what follows, we will present four curves in all figures. They 
correspond to two values of ccp trajectory and presence or absence 
of the charmed sea contribution: 

1) = —2.18 (full line) with corresponding values of param¬ 
eters a q = 0.0007, af = 5, b Ds = 1.6, 5 C = 0.05, 

2) = 0 (dashed line),a 0 = 0.0005, a\ = 5, b Ds = 1.6,5 C = 
0.05, 


3) ap = —2.18 (dashed-dotted line) - the same as in 1) without 
charmed sea contribution S c = 0 and 

4) cp/, = 0 (dotted line) - the same as in 2) but S c = 0. 

The charmed sea suppression factor 5 C = 0.005 was obtained 
earlier in 0. 

All theoretical curves in the model under consideration are 
sums of the directly produced Zl s -meson cross section and the 
contribution from the decay of the corresponding D* —*• D'y res¬ 
onance. 

The experimental data on the xp— dependence of the differ¬ 
ential cross section of LU-meson in Ti~p- interaction at initial 
momenta Pl = 230 GeV [|5j| together with model calculations are 
presented in Fig.l. Note that the normalization of the experimen¬ 
tal data on this figure is arbitrary and here we may compare only 
the shape of the cross section. 

The model calculations for the spectra of the sum of all D s - 
mesons in the reaction n~p —*• D s X at Pl = 350 GeV/c are 
compared with the experimental data at Pl = 350 GeV/c 0 
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in Fig.2. It is necessary to mention that in the irp- collision the 
production cross sections of D~(sc) and Df(cs)~- mesons are equal 
because they do not contain the valence quarks from pion beam. 
The experimental measurements of the D~ Df- mesons total 
cross sections , given in slightly differ each other, although 
they are close within the rather large experimental errors. 

In Fig.3 the experimental data on the total cross section ver¬ 
sus momentum of the initial i r meson in xp > 0 region are com¬ 
pared with the model calculations. As we can see the existing 
experimental data do not allow to draw a definite conclusion on 
momentum dependence or, as it was noted in 0, there is no indi¬ 
cation of momentum dependence which contradicts to PYTHIA 
(see curve on Fig.9c in |7]]) and our calculations. It is also neces¬ 
sary to note that the first point at Pl = 235 GeV/c (NA32) || 
stands for the D+ meson cross section, while two others - (E769) 
at Pl = 250 GeV/c jT(| and (WA92) at Pl = 350 GeV/c 0 were 
given for the sum of Df and D~ mesons. 

The predictions of the model for the inclusive spectra of D+ 
mesons in n~p- collision at Pl = 500 GeV/c and T,~p collisions 
at Pl = 330 GeV/c and 600 GeV/c are given in Figs.4 - 8. 

The xf dependence of leading/nonleading production asym¬ 
metry for Y,~p- collision together with experimental 

data of the WA89 collaboration for P L = 340 GeV/c 
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is 

presented in Fig.9. The predictions for the same asymmetry at 
Pl = 600 /c is shown in Fig. 10. 

We do not present the comparison of our calculations for asym¬ 
metry (D~,D+) in 7 T~p collision with the experimental measure¬ 
ments |T2| because it is equal to zero due to equality of D~ and 
Df mesons cross sections. 

As we can see from the comparison with the experimental 
data, the model under consideration represents the shape of D s - 
rnesons cross section in 7i~p- collision. Slightly better agreement 
is for the case of charmed trajectory intercept equal to = 0. 
Recent experimental measurements of leading/nonleading asym- 
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metry (Dj,D+) in the E _ p- interaction at Pl = 340 jc |7J prefer 
the = —2.18 without charmed sea contribution (curve 3 on 
fig.9). Unlike the 7r induced reactions where, as it was already 
noted, the asymmetry is equal to zero, in the hyperon induced re¬ 
action it should be significant due to the presence of the valence 
s- quark which was testified by WA89 collaboration measure¬ 
ments 0. However in the description of inclusive cross section 


the variant 3) gives the rapidly falling cross section. Taking into 
account the charmed sea contribution leads to decreasing of the 
cross section, especially at large xp, but leads to de decreasing of 
the asymmetry at the same region. Unfortunately existing exper¬ 
imental data due to absence of the cross sections measurements 
at xp > 0.5 do not allow to do the unambiguous choice between 
considered variants. 

In the model under consideration it is possible to improve 
the description of the experimental data [0] by slightly changing 
the value of the cq parameter in the fragmentation functions of 
the s- quark (|3|) and ds- diquark (||). However for the exact 
determination of the quark and diquark fragmentation functions 
into the D s - mesons more precise measurements of all observa¬ 
tions are needed: total and differential cross sections and lead¬ 
ing/nonleading production asymmetry especially at large xp. 
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Figure captions 

Fig.l. Comparison of the QGSM calculations with the ex¬ 
perimental data (NA32) || on differential cross sections of the 
Dt -meson production in 7r p- interaction at Pl = 230 GeV/c. 
The theoretical curves notification: 1) full curve were calculated 
for otf = —2.18, 2) dashed curve stands for — 0 (in 1 and 2 
with the contribution of the charmed sea), 3) dashed - dotted line 
- — —2.18, 4) dotted line - = 0 ( in 3 and 4 the charmed 

sea contribution is absent). 


Fig.2. Comparison of the QGSM calculations with the exper¬ 
imental data (WA92) 0 of the summary spectra of H s -mesons 
production in Ti~p- interaction at 350 GeV/c. The curves are the 
same as in Fig.l. 


Fig.3. Comparison of the QGSM calculations with the experi¬ 
mental data on D s mesons total cross sections in 7i~p- interaction. 
The curves are the same as in Fig.l. 


Fig.4. Predictions for the differential cross sections of the in¬ 
dependence of the TC-meson production in n~p- interaction at 
Pl = 500 GeV/c (E769). The curves are the same as in Fig.l. 

Fig.5. Model prediction for the in-dependence of the D~- 
meson cross section in interaction at 340 GeV/c. The curves 
are the same as in Fig.l. 


Fig.6. Model prediction for the x /^-dependence of the Df- 
meson cross section in interaction at 340 GeV/c. The curves 
are the same as in Fig.l. 
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Fig.7. Model prediction for the xf~ dependence of the D~- 
meson cross section in E“p interaction at 600 GeV/c. The curves 
are the same as in Fig.l. 


Fig.8. Model prediction for the ^-dependence of the -D+- 
meson cross section in interaction at 600 GeV/c. The curves 
are the same as in Fig.l. 


Fig.9.Comparision of the model calculations for the a^-dependence 
of D~/Df asymmetry for the beam at 340 GeV/c with the 
experimental data |ll|. The curves are the same as in Fig.l. 


Fig.10. Predictions for the rcp-dependence of D~/Df asym¬ 
metry for the E _ beam at 600 GeV/c. The curves are the same 
as in Fig.l. 
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